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Abstract. The present paper deals with the homogenization of the heat 
conduction which takes place in a binary three-dimensional medium consisting 
of an ambiental phase having conductivity of unity order and a rectangular 
honeycomb structure formed by a set of thin layers crossing orthogonally and 
periodically. We consider the case when the conductivity of the thin layers is 
in inverse proportion to the vanishing volume of the rectangular honeycomb 
structure. We find the system that governs the asymptotic behaviour of the 
temperature distribution of this binary medium. The dependence with respect 
to the thicknesses of the layers is also emphasized. We use an energetic method 
associated to a natural control-zone of the vanishing domain. 
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1 Introduction 



The study of the lattice-type structures which are characterized by period- 
icity and small thickness of the material is one of the main achievements of 
the homogenization theory. Civil engineering, electrotechnics and the aerospace 
industry are mainly concerned by composite materials, particularly truss struc- 
tures, which have to be treated by techniques of homogenization when direct 
computations fail. 

The foundations of homogenized layered materials were laid down by Murat 
and Tartar in their pioneering work [lOj . This method was still used by [9] in 
the framework of weaker topologies. The case of i?V^-functions and sequences of 
measures is worked out in [6] . Many examples and applications may be found in 
[T] . The difficulties arising from the direct computation of the behaviour of these 
structures are twofold: the great number of cells and the small thickness of the 
material. The periodic distribution classically suggests that the homogenization 
method for perforated domains should be used as in [11] , where the period e 
and the thickness 6 of the material are considered as independent vanishing 
parameters. A classical issue is then to consider the problem of permuting both 
converging processes. This is made in details in As the pratical point of 
view favors geometric consideration, measures provide a more general tool when 
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they are confined to the description of the critical part of a system, as in [5] 
and [3]. This approach was used in |12| . where the homogenization of elastic 
reticulated structures is performed for an anisotropic material surrounded by 
an empty environment. 

Our study is based on the control-zone method which was introduced and 
developed by [1]-IS] and I7]-[3, specific to the binary composition of the system, 
namely the truss and the ambicntal phase. This procedure proved its efficiency 
in the modelling of fine substructures where concentrated material of high con- 
tribution influences the behavior of the global problem in spite of its vanishing 
volume. The asymptotic treatment reveals the apparent paradox between an 
obviously disappearing element and its everlasting action on their environment. 
The coupling between both components manifests itself through the rarefying 
ratio 7e which arises as a criterium for the reduced problem to exist. Unlike the 
critical case studied in [5] , the connections between the layers of the reticulated 
structure annihilate the capacitary term lime^o 7e when it is defined. This was 
also observed in [5] in the absence of connections. Interestingly, the presence 
of the ambiental phase substantially modifies the contribution of the truss in a 
way that cannot be anticipated from [11]- [12], because the sequence of problems 
under consideration behaves singularly with respect to the period of distribution 
and involves a new criterium, namely the capacity of the intersections of the 
crossing layers. 

The present paper is organized as follows: 

Section [2] is devoted to the main notations and to the description of the 
initial problem. We set the functional framework for which the existence and 
uniqueness of the solution can be established. 

In the first part of Section [3] we present two operators that have a localiz- 
ing effect and that we use in order to obtain the specific inequalities related to 
honeycomb structures. We also define the capacitary functions and step approx- 
imation operators which are associated with the test-functions defined in our 
control-zones. They overcome the singular behavior of the energy term when e 
tends to zero. 

Section [J] deals with the homogenization process in the reticulated case in 
the so-called box structure geometry. We obtain the homogenized equation 
which displays explicit effective coefficients. The proof relies on the energetic 
method applied in the control-zone context. No critical thickness of the layers 
appears and the intersections of the layers have no distinct infiuence upon the 
asymptotic behaviour of the temperature. 

The homogenization process in the gridwork case is worked out in Section (5] 
We obtain the homogenized equation and explicit effective coefficients in this 
case also. The proofs are only sketched because the arguments follow the same 
lines as in Section H] 

The significant difference between the effective coefficients obtained in the 
two cases studied here shows how important is the internal geometry of the 
vanishing superconductive material. 
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2 Setting of the problem 



Let ft = , I =] , - [ and n £ N. From now on we use the notations 

£=^-^^, Z. = {fc e Z, \k\<n\ and I^=ek + el, keZ^. (1) 

Obviously, cardZ^ = 1/e and x G / if and only if there exists fc G Z^ such 
/j. 

that X G I^. 

ri 

For any i G {1,2,3}, we consider rl > 0, rl « e, that is — — > when 
£ ^ 0. For any fc G Z^ we define: 

T^ fe = {a; = (xi,X2,X3) G n, \x, - sk\ < rl}, = U^^zTlk (2) 

= uLiTe, = n r/, for l < ^ < j < 3. (3) 

We consider that Q is occupied by two materials with highly different conduc- 
tivities; one is highly conductive and it is concentrated in the vanishing domain 
of the thin layers {\Tg\ — > 0) and the other forms the ambicntal phase which 
has conductivity of unity order. 

Assuming that we are given a source term in fl and some conductivity 
coefficients a,b > 0, we consider the heat conduction problem, that is we study 
the temperature Ue which satisfies in some sense the system 



- aAu^ = A in n\T, (4) 
-^A«, = A in T, (5) 
= on dil, (6) 



together with the natural transmission conditions on the interface dT^ \ dil. 

More precisely, assuming that fe G H^^{il.), the variational formulation of 
our problem is : 

To find Me G H^{Q) such that 

a I \/ue\/v + b-f VueVv ^ {fe,v), VuGi?o(17), (7) 
Jn\T, J 

where we have used the notation 

/• = / ■ for any measurable E C fl (8) 
E \E\ Je 

and (•, •) is the duality product between H~^{ft) and Hl{n). 
Applying the Lax-Milgram theorem wc obtain: 

Proposition 2.1 The variational equation has a unique solutionu^ G H^(fl). 
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The aim of this paper is to describe the asymptotic behaviour of the tem- 
perature Me as £ — > 0, assuming that the source term is weakly convergent: 

3feH-^{n) such that fe^f in H^^{n). (9) 

We have to remark here that the boundedness of imphes immediately: 

Proposition 2.2 (u^)^ is bounded in HQ{fl) and there exists C > 0, indepen- 
dent of e, such that 

|VueP<C. (10) 



3 The control-zone homogenization method 

In order to obtain further results, specific to the thin substructure considered 
here, we have to introduce the following operators: 

Definition 3.1 To any u G H^{Vl) we associate G\{u) G i^(ri) defined by the 
following: 



(••■, A, (11) 



fceZe 



Xi—ek — rl. 



\xi—6k-\-rl 1 



fc G Ze 



(12) 



where 



^k±ri the traces of u on the corresponding faces of ^ . 



These operators have three basic properties, which were already proved in 

[HI: 

Proposition 3.2 For any u G H^{Vl) and i G {1,2,3}, we have 



\Gl{u)~u\' <rl 



di 



dxi 

\Gl{u)\^ = \Gl{u)\l 
du 



dxi 



\Gl{u)-u\n<e 

where \ ■ \n is denoting the norm of L^(fl). 

The first important consequence is 
Theorem 3.3 There exists G > 0, independent of e, such that 



(13) 

(14) 
(15) 



(16) 
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Proof. Let us notice that 

2 2 

/ \u\'<J2i I"P<2E/ ,{\Gl{u)^u\^ + \Gl{u)\^). 
Using (HSll-llll), it yields 

2 

up < 2 ( max^r^) \^u\l + 4^] {\G\{u) - u\l + ) 
and the proof is completed by (|15p and the Poincare-Friedrichs inequality in VL. 



As the techniques of Section 3 [5] can be used to the domain T^'-* (i < j), 
then, according to [2], we have: 

Theorem 3.4 There exists C > 0, independent of e, such that 

\u\' <Cma.Jl,e^\n^,e^\n-\\Vu\l, VueH^{n). (17) 



Finally, we remind an estimation of the same type proved in [5] and which 
is associated to 

T° = nf^iT;. (18) 

Theorem 3.5 There exists C > 0, independent of s, such that 

j- < CmaxJ^l, |Vu|^j, \JuC,Hl{n). (19) 

As the vanishing highly conductive layers have thicknesses of size < 
r\ << £, we find the asymptotic behaviour of by applying a control-zone 
method, that is an energetic method using test-functions associated to an ad- 
equate control-zone: a vanishing set which includes the layers and has much 
larger thicknesses i?*, where rl « R\ « e . 

First, let us introduce the tools of this method. Denoting 

n^{{Re)e, r,«R,«e}, (20) 
then, for any {RV)e € TZ wc define the control-zone of the present problem by 

C, = utiC:, C: = Ufcez.Qfc, *e {1,2,3}, (21) 
where for any fc G and i G {1, 2, 3} we have 

Clk^{x^{xi,X2,X3) en, \x^-ek\<Rl}. (22) 

The test-functions associated to this control-zone are defined by using the 
following capacitary functions wl G W^'°°{fl) {i = 1,2,3), given by 



wl{xi,X2,X3) = < 



K 

1 - if .T = (xi, X2, xz) G Tl 

1-^^ if .T G(q,,\T;,,) for some fcGZ, ^^^^ 

^ if .tgO\Q, 
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and the step approximation operators introduced by 

Definition 3.6 To any Lp G T>(fl) we associate (fil G L°°{n) (i = 1,2,3), as 
follows 

(pl{xi,X2,X3) = (p\^.^=^k{Xt) Ijk (Xi), (24) 

where ek + 2WJ . 

These operators have the foUowing basic properties: 

2X1/2 



|V<|c. < ( ^ ) (25) 



|(^-(^^U^(cj) <i?^|V(^U^(o). (26) 



1 /2 

|V<^^|c. < ( ^) |V^|i=.(n). (27) 



4 The reticulated case 

For any i G {1, 2, 3} let us denote m,; > as the hmit of 

"^. = limS- (28) 

\Te\ 

Obviously, we have 

nil + ^2 + iT^z = 1- 
In this section, we consider the case when 

>0, ViG {1,2,3}, (29) 

that is the case when all the three parameters r* have the same order of magni- 
tude with respect to e. This geometry is called sometimes as the box-structure 
case. 

We can present now the preliminary convergence results: 

Proposition 4.1 There exists u G H^(Q,) such that, on some subsequence, 
there hold 

u^-^u in H^{n) (30) 

Gl{u,)^u in L\n), ViG {1,2,3} (31) 

/ u,v^ ( uv, VvEH^in). (32) 
Moreover, for any i G {1, 2, 3} we have: 

A.lj"""!^''' ^''^^o^")' ^J'^il'^.S}, jV*- (33) 
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Proof. The weak convergence pop follows from Proposition 12.21 As i?o(r2) is 
compactly embedded in L^(fJ), the strong convergence pTjl is obtained by using 

dUD and dsni). 

In order to prove (j32p we remark that 

uev^^-f-l u^v- .-^- + E^f (34) 

n J TP i<f^<3 J T- ^ \Te\J Ti 

Using ([TU)) . (fT7|) and (fTO|) . we prove that the first two terms of ([M)) are converging 
to zero, as follows: 



\Te\Jl 

f u^v 



< C- 



■ max 1 



ITU I 



,3/2 
J72 



1 



max 1, eln^— , eln^- 



|Vw|o -^0, 
1 



n 



Taking into account ([28|) . the proof of ([32|) is completed by applying to Propo- 
sition 2.15 of where we have proved that 



,v~^ [ uv, Vv e H^{n). 
Jn 



(35) 



Next, let (f e T>{n) and i G {1,2,3}. Denoting hy v — {vi, 1^2, 1^3) the outward 
normal to dT^ we obviously have (piyj = on dT^ 7^ *) and hence: 



/du^ f dtp 



Using ((35)) . from the previous relation it follows 



dx 



-ip 



dp 
dxi 



du 



n dx. 



■p. 



The proof of ([55]) is completed by continuity, using (fTU]) and ([T5)) . ■ 

Now we can present our main result. 

Theorem 4.2 (we)e *s weakly convergent in HQ{fl). Its limit, u G H^{Q,), is 
the only solution of the equation 



f in fl. 



(36) 



Proof. For any i S {1, 2, 3}, let (Rj)e G TZ and p G P(rj). Using the definitions 
(I23l)-(l2ll, we denote 



(37) 



7 



We notice that 



iplix)wl{x) = < 



1 - 



\xi - ek\ 



v\xi=ek[xi), X ^Tl f^ for some fc £ 



(38) 



0, a; e 17\Ci 

and hence v^{}p) € WQ^'°°(ri). Then we set v = v^[ip) in ([7]) and it follows 



X] ( ^ m 



(39) 



As \Ce\ -^0 it follows that 



Rl 



(V(^)lj^\C, ^ strongly in L^{VL) 



(40) 



and using ([50]) we obtain the convergence of the first left-hand side term of ([55]) : 



(41) 



The second left-hand side term of ([M]) is converging to zero. Indeed, as (ug)^ is 
bounded in Hl{^), we have 



a I VUeVWe((/3) 



(42) 



Moreover, by a straightforward computation we find that 

3 

|V^;,(^)|c. < C|V<^|c. +CY,{WI - <p|l~(c-)|V<|c. + |V^.(pi|cJ (43) 



and the assertion is proved by using ((25l) - ((2^ 



|Vz;,(^)k <C|V^U^(o)E(? 



0. 



(44) 



For the convergence of the third left-hand side term of ([55)1 we notice that 

1 



1 



/ '^XiUe{x)VxMxi=ek{xi) dx 



(45) 
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Using ([TU)) and the smoothness of Lp we have 



1^ X! \ V;F,'Ue(x) (V^.cp|x.=efe(a;i) - Vaj,(p(a;)) da; 

A; G Z £ s,k 



< 



and thus we replace (|45|) by 



hm + Vu^Vpl ~ hm 



m 



el JT. 



Next, recahing and applying ([55]) we obtain 



3 

= 65^(1 -mO / 



9u dp 

Finally, noticing that ((25)) and ([26l) imply 

|(/„(^^-<^K}| <C|<|cjV^U^(o)+C|V<|cJ^^-<pk..(Cj) < 

0, 



(46) 





Rl 


l/2\ 




e 





we obtain the convergence of the right-hand side term of ([5^ . that is: 

(/„ v,{p)) =^2(1-^) ifcP) + - - 3 (/, p). (47) 

Resuming, we can say that if we pass ([39]) to the limit, then we obtain: 
3 



3a [ VuV(/7 + 6V(l -m,) [ 
Jn Jn 



du dip 
Q dxi dxi 



3{f,p>), ypev{n), (48) 



which is. by continuity, the variational formulation of ([36]). The convergence of 
the entire sequence (ue)e follows from the unicity of the solution of 



Remark 4.3 In spite of the vanishing volume of the rectangular honeycomb 
structure the homogenized behaviour is anisotropic, except the case when mi = 
1 

m2 = ma = -. 

Remark 4.4 The reticulated domain occupied by all the intersections of the 
layers has no distinct influence upon the asymptotic distribution of the temper- 
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ature. For i < j and k G this follows from 



\Te\ JTl- 



<C|V.|...„,|J^/^JV..P)'(M)'< ^^^^ 



<cm(|)' ^ 0, 

where we have used hlOi). 



5 The gridwork case 

Using the notations and definitions of tlie previous sections, the gridwork case 
corresponds to the situation when the horizontal layers are missing and the 
other two layers have the same order of magnitude with respect to e. It follows 
that 

= Tl U (50) 

and there exist 

\rpi\ 

rrii = lim — ^ (i — 1, 2), such that mi + m2 = 1. (51) 
£^0 I Tel 

For consistency with the honeycomb case, we also define here 

mg = 0. (52) 

Obviously, for i G {1, 2} th properties of the corresponding operators defined 
in Section [3] are still valid. Then we prove results similar to those of Proposi- 
tion [1711 with only one mentionable distinction: 

For any i G {1, 2} we have ([55)) . 

Consequently, the homogenization result in this case is the following: 

Theorem 5.1 {u^)^ is weakly convergent in iJQ(f2). Its limit, u G H^(Q.), is 
the only solution of the equation: 

-Ef«+|(l--0)|^ = / » ^- (53) 



1=1 



Proof. The proof is completely analogous to that of Theorem 14.21 with the 

difference that the control-zones are oriented in only two directions. That is, 
with the test-function 

= E f f 1 - ^ + ~ (54) 



1=1 

we get the same convergence as in (j46p 

2 / i X „ 2 



.=iV ReJJn Jn ^gg^ 



3 

65^(1 -m,) 

i=l 



du dip 
dxi dxi 
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where we have taken in account the convention ([5^ also. From here the proof 
can be completed identically. ■ 



Remark 5.2 Assuming that the quantity \Ts\ is the same in the both cases that 
we considered, the significant difference between the equations i36\) and JJ^p 
shows how important is the internal geometry of the vanishing superconductive 
material. 
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